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ABSTRACT 

It was recently discovered that the temperature in the surface layer of externally heated optically 
thick gray dust clouds increases with the optical depth for some distance from the surface, as opposed 
to the normal decrease in temperature with distance in the rest of the cloud. This temperature 
inversion is a result of efficient absorption of diffuse flux from the cloud interior by the surface dust 
exposed to the external radiation. A micron or bigger size grains experience this effect when the 
external flux is of stellar spectrum. We explore what happens to the effect when dust is a mixture 
of grain sizes (multigrain) . Two possible boundary conditions are considered: i) a constant external 
flux without constrains on the dust temperature, and ii) the maximum dust temperature set to the 
sublimation temperature. We find that the first condition allows small grains to completely suppress 
the temperature inversion of big grains if the overall opacity is dominated by small grains. The second 
condition enables big grains to maintain the inversion even when they are a minor contributor to the 
opacity. In reality, the choice of boundary condition depends on the dust dynamics. When applied 
to the physics of protoplanetary disks, the temperature inversion leads to a previously unrecognized 
disk structure where optically thin dust can exist inside the dust destruction radius of an optically 
thick disk. We conclude that the transition between the dusty disk and the gaseous inner clearing is 
not a sharp edge, but rather a large optically thin region. 

Subject headings: accretion, accretion disks — circumstellar matter — dust, extinction — stars: 
pre-main-sequence 



1. INTRODUCTION 

Dust is one of the principal components of interstellar 
and circumstellar matter. It serves as a very efficient 
absorber of starlight, which is dominated by visual and 
ultraviolet photons, and reemiter of this absorbed energy 
into infrared (IR), thereby modifying the entire spectral 
energy distribution. Hence, radiative transfer in dusty 
environments is an inherent part of the study of star 
formation, protoplanetary disk evolution, dusty winds 
from AGB stars, circumnuclear environment of AGNs, 
etc. 

In general, radiative transfer is extremely complicated 
because the dust in these environments always comes as 
a mixture of various grain sizes (multigrain 1 ), shapes 
and chemical compositions. It is, therefore, common 
to employ certain approximations that make the prob- 
lem computationally manageable. For example, dust 
grains are usually approximated with spheres of simi- 
lar chemical composition, so that the Mie theory can 
be easily used for calculating dust cross sections. An- 
other most commonly used approximation is replacing a 
mixture of dust grain sizes by an equivalent (i.e. syn- 
thetic or average) single grain size. Numerical calcula- 
tions have shown that this approximation does not pro- 
duce a significant change in the spectral energy distri- 
bution (lEfstathiou fc Rowan- Robinsonlll994t IWoldl2003l : 
ICarciofi et al.ll2004l ). It is assumed that this is particu- 
larly appropriate for optically thick dust clouds. 

It has been discovered recently that externally heated 

1 Typically, the term multigrain also includes all other dust grain 
properties, like the grain shape and chemistry. But, for simplicity, 
we use this term only to describe grain size effects. 



optically thick clouds made of large (> 1/xm) single 
size dust grains produce the effect of temperature in- 
version where the maximum temperature is within the 
dust cloud 2 , at the visual optical depth of Ty ~1, in- 
stead of on the very surface exposed to the stellar heating 
(jlsella fc Nattdl200arVffikovic et al.f 2006). The cause of 
this inversion within the surface layer is its ability to ef- 
ficiently absorb the diffuse IR radiation originating from 
the cloud's interior - a process similar to the "greenhouse 
effect." 

It is not known, however, if multigrain dust would pro- 
duce the same temperature inversion effect. Here we em- 
ploy analytical multigrain radiative transfer to study the 
surface of optically thick clouds. Two possible types of 
boundary conditions are explored: i) a constant exter- 
nal flux heating the cloud, with no limits on the dust 
temperature, and ii) a fixed maximum dust temperature 
corresponding to dust sublimation. In we describe the 
analytic method and apply it on single size dust. In the 
next section [J3]we apply the method on multigrain dust. 
After that in 21 we explore a possibility of surface ther- 
mal cooling in transverse direction. Some aspects of our 
result are discussed in $5j with our conclusion in $6l 

2. OPTICALLY THICK SINGLE SIZE GRAIN DUST CLOUD 

In order to understand effects of multigrain dust on the 
temperature structure of optically thick clouds, we first 
have to take a look at the single size dust grain clouds. 

2 Another type of temperature inversion has been recog- 
nized in protoplanetary disks, where additio nal viscous heating 
can increase the disk interior temperature (Calvct et afj 119911 ; 
IMalbet fc Bcrtout 199lJ). In contrast, the temperature inversion 
discussed here is a pure radiative transfer effect and does not re- 
quire any additional assumption (like disk viscosity) to operate. 
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Here we reiterate analytic techniques described in th e lit- 
erature (jlsella fc Nattal 120051 : IVinkovic et al.ll2006| ) and 
then expand our approach to multigrain dust in the fol- 
lowing sections. 

Consider dust grains at three different locations (Fig- 
ure [1]): on the very surface (point Po), at optical depth 
Ty (point Pi) and at optical depth 2ry (point P 2 ). The 
optical depth Ty is defined at the peak wavelength of the 
bolometric temperature of external flux F in illuminating 
the cloud. We will be interested in cases where dust can 
reach sublimation temperatures (between ^1,000K and 
2,000K for interstellar dust) and where the external flux 
source is a star-like object peaking its emission at submi- 
cron wavelengths. For the purpose of this paper we use 
the optical depth Ty in visual, but one can adjust it to 
the required wavelength of interest. 

The surface layer between points Po and Pi emits the 
thermal flux F sur \. We assume that it emits equally 
on both sides (toward the left and right in Figure Q]) . 
This approximation is valid if the layer is optically thin 
at the peak wavelength of its IR emission. If we de- 
fine q = cry/criR, where ay and am are dust absorption 
cross sections in visual and IR, then the IR optical depth 
requirement is 

rm = — < 1. (1) 

The same requirement holds for the second surface layer 
between points Pi and P2 and its thermal flux F sur 2- 
Each of these two layers attenuates the external flux 
by exp(-Ty) and the thermal fluxes by exp(— tir) — 
exp{-T V /q). 

Since the cloud is optically thick to its own radiation, 
no net flux can go through the cloud. The total flux 
coming from the left at any point in the cloud is equal 
to the total flux coming from the right. This yields flux 
balance equations at points Pq, Pi and P2 
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where F out is the thermal flux coming out of the cloud 
interior at P 2 . 

We want the external flux Fi n at point P2 to be atten- 
uated enough to make diffuse flux the dominant source 
of heating. This requirement means that 2ry > 1, which 
in combination with equation [T] gives the allowed range 
for Ty 

0-5 < t v < q. (5) 
From equations [2HH we further get 
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A dust grain at Po absorbs ~ ayFi n of the external di- 
rectional flux. It also absorbs ~ 2amFjn of any infrared 
diffuse flux Fir, where the factor 2 accounts for absorp- 
tion from 2-7T sr. The grain emits at its temperature T 
into 47r sr, so that the energy balance is 

a v F m + 2a IR (F surl + F sur2 + F out er 2Tv ' q ) = 




Fig. 1. — Sketch of the surface of an optically thick single size 
j;rain dust cloud illuminated from outside. See [[2]for details. 
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where asB is the Stefan-Boltzmann constant. Using the 
flux balance in equation [2] and q = ay /am, we get 
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Similarly, we can write the energy balance for a dust 
grain at Pi and P 2 

2F sur2 + 2F out e~ Tv / q = Aa SB T* 

(10) 
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We approximate the interior flux as 
Fout ~ vsbTo- 



(12) 

This is a good approximation of the interior for gray dust 
and an overestimate for non-gray dust, where tempera- 
ture decreases with optical depth. Now we can continue 
deriving temperatures in two possible ways. 

METHOD 1: Combining equations EE |9l [TUj \U\ and 
QJ yields 
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B = qe- 2Tv + 2(1 - e - Tv - Tv,q ) 
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The upper panel of Figure [2] shows this temperature 
profile for Ty — 1. The effect of temperature inversion 
T 2 > To is present when q < qu m u ~ 1.4. Only grains 
larger than about one micron can have such a small value 
of q, which makes them almost "gray" in the near IR. 
This inversion does not appear in non-gray dust (q > 
qiimit) where the temperature decreases monotonically 
with distance from the cloud surface. 

METHOD 2: We can use equations HH] to express F out 
as a function of Fi n and then use equations [9] and [12] to 
obtain 



T 9 4 = 



1 



(1 



'/q\ 



1 -|- e -rv/q 



(17) 



Similarly, we can express F aur \ and F sur 2 as a function 
of Fi n and then use equation [TU] to obtain 
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where q\ = o~v.il cira- This shows that dust grains of dif- 
ferent sizes have different temperatures on the surface of 
an optically thick cloud. The lowest temperature among 
the grain sizes is acquired by the largest grains because 
they have the smallest qi. 

On the other hand, the contribution of the external flux 
is negligible Fi n exp(—Tv) — > when ry > 1, therefore 
F SU r — Fouti which yields the same temperature for all 
grain sizes 

T?{T V ^l)=F out /o~ SB . (21) 

A general relationship between grain temperatures is 
derived by subtracting equation [19] for a grain i from the 
same equation for a grain j and then use equation [^0] to 
obtain 



-e- TV Tf{Q) (22) 



tHt v )=T 1 Ht v ) + 



Since q scales inversely with the grain size, equation 1221 
shows that smaller grains always have a higher temper- 
ature than bigger grains at any point in the cloud, with 
the limit Tj ~ 7j when exp(-Ty) <C 1. 

Optical depth is now a cumulative contribution of all 
grain sizes in the mix. 



JV 



(23) 



visual optical depth 

Fig. 2. — Surface temperature of an optically thick single size 
grain dust cloud for various visual to infrared opacity ratios q. 
The upper panel shows the solution based on equations 1131 and 
II II (method 1). The lower panel is based on equations 1171 and 1181 
(method 2). Optical depth step is Ty = 1. Both solutions show 
temperature inversion (increase of temperature with optical depth) 
for small q. 



The lower panel of Figure [2] shows this result for ry = 1 . 
It differs slightly from previous solution in Method 1, 
but gives qualitatively the same result. The temperature 
inversion now exists for q < qu m it ~ 1.5. Henceforth, we 
use Method 1 for further analysis. 

3. OPTICALLY THICK MULTIGRAIN CLOUD 

3.1. Constant external flux as the boundary condition 

We consider now an externally heated optically thick 
dust cloud made of N grain sizes, under the condition 
of constant external flux and no limits on the dust tem- 
perature. A dust grain of the i th size at a point of ry 
optical distance from the cloud surface is heated by the 
attenuated stellar flux Fi n exp(—Ty) , by the local diffuse 
flux F sur coming from the direction of the surface, and 
by the local diffuse flux F out coming out of the cloud in- 
terior. The thermodynamic equilibrium of the dust grain 
gives 

o~v,i Fi n e~ TV + 2a- JR j (F sur + F out ) = 4:0- IRti o- SB T t (ry). 

(19) 

Factors 2 and 4 account for absorption from 2ir sr and 
emission into 47r sr. 

On the very surface of the cloud ry = and F sur = 0, 
therefore Fi n — F out , which yields 
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where nt is the number density of the i th grain size. If 
we scale optical depth relative to ry then 
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where Ta,« is the relative contribution of the i th grain to 
the dust opacity at wavelength A 
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Consider now a model similar to figure [T] except that 
the dust is multigrain. For simplicity and clarity of the 
following analysis, we work with only two grain sizes: a 
"big" grain a with q a ~ 1 and a "small" grain f3 with 
5/3 > Qa ■ The infrared optical depth step is 

TIR = T V 1 ■ (26) 

V la qp J 

According to equation l22l the temperature of small grains 
at optical depth fcry (point P^) is 
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The temperature of big grains at point Po is (equation 

"IB 



F 
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4ctsb 

while for the other temperatures we need the flux balance 
at points Po, Pi and P2. Since the balance is the same 
as in equations [2][4j except that the infrared step ry jq 
is replaced by tir, we use the procedure described in 321 
and derive 

Tt, a = -^TL + ^L (29) 
q a + 2 la sb 
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Fig. 3. — Temperature of two types of grains in the surface of an 
optically thic k multi grain dust cloud: larger grains a {thick l ines , 
equations !36l and l37H and smaller grains /3 (thin lines, equation l27l l. 
Line styles correspond to different ratios q a of visual to infrared 
opacity of big grains. The same ratio for small grains is fixed to 
9/3 = Qa + 1. Panels show results for two relative contributions of 
big grains to the visual opacity: Tv a = 0.1 in the upper panel and 
Tva = 0.9 in the lower panel. Optical depth step is ry = 1. Notice 
how the temperature inversion in big grains (small q a ) disappears 
when small grains dominate the overall opacity (upper panel). The 
inversion exists only if bigger grains dominate (lower panel) . 
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Fig. 4. — Same as Fig. [3] but for qp = q a + 10. 
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The interior flux F out is now a cumulative contribution 
of all grain sizes according to their relative contribution 
to the dust opacity. In our two-size example 

F ou t jo SB — f/fl,ce^2,a + IR,fi^2,p- (33) 

Combined with equation [23 gives 



Font jo SB = T 2 
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Putting together equations [29l [30] and [34] yields the 
solution 

rpA _ Ag + Bg + Cg ^4 
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The resulting temperature is plotted in Figures [3] and 
[U The upper panel in each figure shows the result when 
small grains dominate the opacity (Ty Q = 0.1). The 
temperature inversion in big grains is suppressed because 
the local diffuse flux is dictated by small grains. If big 
grains dominate the opacity (lower panels, Ty, a = 0.9) 
then the temperature inversion is preserved. 

3.2. Sublimation temperature as the boundary condition 

Now we consider a dust cloud hot enough on its illu- 
minated surface to sublimate dust grains warmer than 
the sublimation temperature T su {,. The flux entering 
the cloud is adjustable to accommodate any temperature 
boundary condition. From equation 1201 we see that small 
grains are the first to be removed from the immediate 
surface. If the external flux is high enough then the im- 
mediate surface is populated only by % ~ 1 grains ("big 
grains" ) . All other grains ( "small grains" ) would survive 
somewhere within the cloud, at a distance where the lo- 
cal flux is reddened enough by big grains to be absorbed 
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Fig. 5. — Sketch of the surface of an o ptica lly thick multigrain 
dust cloud illuminated from outside. See i|3.2l for details. 
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Fig. 6. — Same as Figure[3] except that small grains are removed 
from the immediate surface by sublimation (Figure [5). The tem- 
perature of larger grains is calculated from equations I43I and I44I 
Despite being mixed with smaller grains, larger grains show tem- 
perature inversion even when smaller grains dominate the opacity. 



less efficiently. 

We again apply our two-size dust model, except that 
the surface layer of visual optical depth ry is occupied 
only by big grains (see Figure [5]). It is too hot for small 
grains to survive within this layer. Both grains exist at 
optical distances larger than ry from the surface. Fol- 
lowing the same procedure as in !T2J we write the flux 
balance at points Pq, Pi and P2 (see Figure [5]) 
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where the IR optical depth between points Pi and P2 is 
given in equation 1261 
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Fig. 7. — Same as Figure[6] but for qp = q a + 10. 



Small grains do not exist now at point P , but their 
temperature at other points is still described by equation 
1271 Deriving big grain temperatures at Pi and P2 is now 
a straightforward procedure already described in $2] and 

go 
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The resulting temperature is plotted in Figures [5] and [JJ 
Two important results are deduced for big grains from 
this solution: i) big grains maintain the temperature in- 
version and ii) if small grains set the maximum temper- 
ature limit then big grains, which are always colder than 
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small grains (see eauation l22[) . cannot reach the sublima- 
tion temperature. 

However, since we do not put limits on the external 
flux, the solution that allows the maximum possible ex- 
ternal flux is the one that also maximizes F out . From 
equation [33] we see that the maximum F out = osBT su b 
is achieved when T2. a = F2.fi = F su b. Such a solution is 
not possible unless Tjfl ;| a = 0. Therefore, dust tempera- 
tures are maximized when all small grains sublimate away 
from the surface region and exist only at optical depths 
of exp(-Ty) « 1. In other words, the external flux is 
maximized when the cloud surface "belongs" exclusively 
to big grains only. 

4. TRANSVERSE COOLING OF THE CLOUD SURFACE 

Solutions presented so far assume no time variability 
in any of the model parameters. The external flux is ad- 
justed by hand to the value that maximizes the external 
flux and keeps dust temperatures below sublimation. In 
reality, however, this is a dynamical process where the 
equilibrium is established by dust moving around and 
sublimating whenever its temperature exceeds the subli- 
mation point. 

The existence of temperature inversion is difficult to 
understand under such dynamical conditions. Since the 
very surface of the cloud is below the sublimation point, 
its dust can survive closer to the external energy source 
than the dust within the cloud. Therefore, the distance 
of the cloud from the source is not defined by the very 
surface, but rather by the dust at the peak temperature 
within the cloud. Imagine now that the whole cloud is 
moving closer to the source. From the cloud point of ref- 
erence, this dynamical process is equivalent to increasing 
the external flux. According to solutions in [J2] and Sj3] 
the peak temperature within the cloud will exceed sub- 
limation at a certain distance from the source and the 
dust will start to sublimate. This point is the distance 
at which the dust cloud seemingly stops. However, there 
is nothing to stop the very surface layer of the cloud to 
move even closer than the rest of the cloud because its 
temperature is below sublimation. 

Notice that we can not resolve this issue in the approx- 
imation of an infinite dusty slab because the transverse 
optical depth (parallel to the slab surface) is always infi- 
nite. No flux can escape the slab in the transverse direc- 
tion. Hence, the radiative transfer in a slab does not de- 
pend on spatial scale and only the optical depth matters. 
The spatial extension of the surface layer is irrelevant and 
solutions from SJ3] and are applicable irrespectively of 
the dust dynamics on the spatial scale. 

In reality, however, the cloud is finite and dust subli- 
mation can eventually make the cloud optically thin in 
the transverse direction. This optical depth gap enables 
thermal radiation to flow out in the transverse direction 
and thus provides a channel for thermal cooling. Under 
such conditions the surface layer cam move closer and 
closer to the energy source by simultaneously expanding 
the size of the transverse optical depth gap, which then 
increases the amount of escaping thermal flux. Since the 
very surface dust is getting closer to the energy source, its 
temperature increases and eventually it reaches the dust 
sublimation. At that moment the radiative and dynam- 
ical equilibria are established and the whole dust cloud 
seemingly stops. 




Fig. 8. — Sketch of the surface of an optically thick dust cloud 
illuminated from outside and cooled through the surface layer in 
the transverse direction. See S|4] for details. 
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Fig. 9. — Surface temperature of an optically thick single size 
grain dust cloud cooled transversely through the surface layer 
(equations I47I and |48> . The amount of transverse cooling is ad- 
just ed to the maximum temperature T(0) at the very surface (equa- 
tion |49)l. Line styles correspond to different ratios q of visual to 
infrared opacity. Optical depth step is ry = 1. Increase in T(0), 
combined with the temperature inversion (see also Fig. |2j, creates 
a local temperature minimum within the surface layer. 



The resulting dust cloud has a large surface zone where 
the thermal flux can transversely escape the cloud. Here 
we make an attempt to describe this on a quantitative 
level through a model shown in Figure [8] It is similar to 
the infinite slab model described in [21 except that now 
we introduce the flux F ex n exiting the cloud transversely 
through the surface layer. 

The flux balance equations [3] and [4] remain unchanged, 
while equation [5] is changed to 

P. — p .. J- P 1 4. P n p—rv/q 4. p , p -2rv/? I AX.\ 
± in ± exit \ ± surl i ± surZ c i ± out ^ ■ V u / 

By following the same procedure as in £j2j we obtain tem- 
peratures Tj! at points 
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(47) 

F exit (48) 



where A and B are given in equations [15] and 1161 The 

values in parentheses are always positive. 
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We compare this with the original solution in equations 
[H [13] and [TU Notice that keeping the interior tempera- 
ture to be the same as before (T!^ =T<z) requires a larger 
surface temperature (Tq > T ). If the surface dust of Tq 
is dynamically driven toward the external energy source, 
as we discussed above, then it will eventually reach the 
sublimation temperature T su b. This dynamical process 
is accompanied by dust sublimation in the cloud inte- 
rior. Sublimation creates a transverse optical depth gap 
and maintains T!{ = T su i,. The established equilibrium 
condition T " = T 2 " = T sub yields 

exlt (q + 2)/(l + e-v/i) - B sub ' { ' 

This flux is positive for values of q that produce temper- 
ature inversion (q < qu m it)- It originates from the diffuse 
radiation trapped within the cloud. Figure [H] shows tem- 
perature profiles for q < qu m it based on F ex n from equa- 
tion 1491 The surface layer now has a local temperature 
minimum at T". 

5. DISCUSSION 

According to equation l2Tj the temperature of dust op- 
tically deep inside an optically thick cloud depends solely 
on the local diffuse flux. Since this can give a misleading 
impression that temperature effects on the cloud surface 
have no influence on the cloud interior, it is worth ex- 
plaining why the surface solution is so important for the 
overall solution. 

At optical depths ry > 2 the flux is exclusively diffuse 
and the dust temperature is dictated by equation 1211 
Notice, however, that even though we set the limit of zero 
net flux flowing through the cloud, there is no limit on the 
absolute value of the local diffuse flux. Arbitrarily large 
but equal diffuse fluxes can flow in opposite directions 
producing the net zero flux. It is the radiative transfer 
in the surface layer that sets the absolute scale. The 
cloud surface of ry < 2 converts the directional external 
flux of stellar spectral shape into the diffuse thermal flux. 
Hence, this conversion determines the value of diffuse flux 
at Ty ~ 2, which then propagates into the cloud interior 
of t v > 2. 

A more rigorous analytic approach shows that the net 
flux is not exactly zero. A very small flux, in comparison 
with the external flux, goes through the cloud and creates 
a temperature gr adient. This yi elds the well known gray 
opacity solution (Mihalas 1978) where T a (t) = ar + b. 
While the constant a depends only on the net flux (a = 
in the case of zero net flux), the constant b is determined 
by the boundary condition on the cloud surface. Hence, 
this becomes a reiteration of the role of the surface layer 
described above. 

The possibility of the temperature inversion effe ct de- 
scribe d in this paper was discussed previously by IWolJ 
(2003). He noticed in his numerical calculations that 
larger grains can have a higher temperature than smaller 
grains in the surface of circumstellar dusty disks. He cor- 
rectly attributed this effect to the more efficient heating 
of large grains by the IR radiation from the disk interior, 
but did n ot analyze it furth er. This is the same effect 
noticed by iDullemo nd (2002) in his numerical models of 
circumstellar disks with gray dust o paciti es. However, it 
was n ot until llsella fc Nattal (|2005l ) and IVinkovic et al.1 
(|2006f ) described this effect analytically in more detail 



that its importance to the overall disk structure was fi- 
nally recognized. 

The main driver for the detailed description of the tem- 
perature inversion effect came from the need for a bet- 
ter understanding of inner regions of dusty disks around 
young pre-main-sequence stars. The advancements of 
near infrared interferometry enabled direct imaging of 
these inner disk regions and resulted in the discov- 
ery of inner disk holes produc ed by dust sublimation 
(M onnier fe Millan-Gabetl 12002] ). However, developing a 
self-consistent model that would incorporate the spectral 
and interferometric data proved to be a difficult problem. 
Two competing models are proposed. In one the data are 
explained by a disk that has a large vertical expansion 
(puffing) at the inner dust sublimation edge due to the di- 
rect s tellar heating of the disk interior (|Dullemond et all 
120011 ). In the other model the inner disk is surrounded by 
an optically thin dusty outflow, without a need for special 
disto rtions to the vertical disk structure (jVinkovic et al.l 
12001 . 

Although current observations cannot distinguish be- 
tween these two models, theory gives some limits on the 
dust properties in the former model. The inner disk 
edge has to be popu l ated by big dust gra ins (> 1/zm) 
(jlsella fc Nattal 120051 : IVinkovic et all 12006ft in order to 
produce a large and bright vertical disk puffing needed 
to fit the data. The hallmark of this radiative trans- 
fer solution is the dust temperature inversion , which we 
also demonstrate in $21 IVinkovic et al.l ((2006) argue that 
purely big grains are unrealistic because dust always 
comes as a mixture of grain sizes (especially in dusty 
disks where dust collisions should consta ntly keep small 
grains in the mix). On the other hand, llsella fc Nattal 
(|2005h point out that smaller grains should be sublimated 
away from the very surface of the inner disk edge, but 
cannot prove that this process is efficient enough to pre- 
serve dust temperature inversion. 

Our analytical analysis in proves that multigrain 
radiative transfer solutions keep small grains away from 
the immediate surface of optically thick disks and pre- 
serve temperature inversion. This is achieved under the 
boundary condition of maximum possible temperature 
reached by all dust grains. Then the surface becomes too 
hot for small grains to survive, which leaves it populated 
only by big grains. By choosing this favorable boundary 
condition, we made an important presumption about the 
dust dynamics. Dust has to be dynamically transported 
to the distances where all grains start to sublimate. This 
can be naturally occurring in optically thick disks due 
to dust and gas accretion. However, in JJUwe discovered 
that big grains can survive closer to the star than the 
inner edge of optically thick disk. The only requirement 
is that the disk becomes vertically optically thin at these 
close distances. 

Therefore, the dust sublimation zone is not a simple 
sharp step-like transition, but rather a large zone (figure 
I10p . We estimate its size by cons idering the inner radi us 
of optically thick and thin disks (|Vinkovic et al.l [2006) 

(50) 

where T*, i?* and are the stellar temperature, radius 
and luminosity, respectively, and \& is the the correction 
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Fig. 10. — Dust dynamics in protoplanetary disks brings dust 
closer to the star than the inner edge ol optically thick dusty disk. 
According to the radiative transfer solution in jj4) such a disk struc- 
ture is possible if this zone is populated only by big grains and it 
is vertically optically thin. These conditions enable disk thermal 
cooling and dust survival (see also figure [Sj . The resulting sub- 
limation zone spans from t he o ptically thin to the optically thick 
inner disk radius (equation 1501 ■ Its physical size is ~ 0.2AU in 
Herbig Ae stars and ~ 0.03AU in T Tau stars. 



for diffuse heating from the disk edge interior. Optically 
thick disks have \t = 2, while optically thin disks have 
\& ~ 1.2. The sublimation zone exists between these two 
extremes, which translates to ~ 0.2AU for Herbig Ae 
stars (~ 5OL ) or - 0.03AU in T Tau stars (~ L Q ). 
This is a considerable size detectable by interferomet- 
ric imaging and has to be addressed in future studies. 
Also, its gas is enriched by metals coming from the sub- 
limated dust, which makes this zone a perfect place for 
dust growth. Notice that, in addition to sublimation, 
grain growth is another way of making the disk optically 
thinner. Hence, the evolutionary role of such a large sub- 
limation zone has to be studied further in more detail. 

Another major problem with the dynamics of big 
grains is that they tend to settle toward the midplanc, 
which suppresses vertical expansion of the inner disk 
edge and leads to the failure of the disk puffing m odel 
(|Dullemond fc Dominik|[200l IVinkovic et al.ll2006D . On 
the other hand, the concurrent model incorporating a 
dusty outflow lacks a convincing physical process respon- 
sible for the formation of a dusty wind. The sublimation 
zone may play an important role in dusty wind processes, 
since its small optical depth and small distance from the 
star should also result in gas properties that are more sus- 
ceptible to non-gravitational forces capable of launching 
a dusty wind (such as magnetic fields) than the rest of 
the dusty disk. 

6. CONCLUSION 

We analyzed the temperature structure of externally 
heated optically thick dust clouds. We focused on the re- 
cently discovered effect of temperature inversion within 
the optically thin s urface of a cloud populated by big (> 
l^tm) dust grains (Ise lla &: Nattal 120051 : IVinkovic et al.l 
2006). The effect is produced by reprocessing the exter- 
nal radiation and not by any additional energy source. 



The inversion manifests itself as a temperature increase 
with optical depth, before it starts to decrease once the 
external directional radiation is completely transformed 
into the diffuse thermal flux. Since small grains do not 
show this effect, the open question was whether small 
grains would manage to suppress this effect when mixed 
with big grains. 

We show analytically that small grains remove the tem- 
perature inversion of big grains if the overall opacity 
is dominated by small grains. However, this does not 
happen in situations where the cloud is close enough to 
the external energy source for dust to start sublimat- 
ing. Small grains acquire a higher temperature than big- 
ger grains and sublimate away from the immediate cloud 
surface. The exact grain size composition of the surface 
depends on the amount of external flux because different 
grain sizes sublimate at different distances from the sur- 
face. These distances are smaller than in optically thin 
clouds because bigger grains provide shielding to smaller 
grains from direct external radiation. We show that the 
temperature inversion is always preserved if small grains 
are removed from only Ty ~ 1 of the immediate surface. 

If the boundary condition requires all dust sizes to 
maximize their temperature, then all small grains are 
removed from the surface layer. They can exist only 
within the cloud interior exp(— ry) <C 1 where the exter- 
nal radiation is completely absorbed. Such a condition 
is expected in protoplanetary disks where dust accretion 
moves dust toward the star. The inner disk radius is then 
defined by the largest grains, no matter what the overall 
grain size composition, because the largest grains survive 
the closest to the star and dictate the surface radiative 
transfer. 

A new problem arises in that case. Since the tem- 
perature inversion keeps the very surface of the cloud 
below the sublimation temperature, its dust can move 
even closer to the star. We show that this creates an op- 
tically thin dusty zone inside the dust destruction radius 
of an optically thick disk (figure [10]). Only big grains 
can survive in this zone. We estimate that its size is 
large enough to be detected by near IR interferometry. 
It consists of big grains and gas enriched by metals from 
sublimated dust, hence favorable for grain growth. This 
shows that the geometry and structure of inner disks can- 
not be determined by simple ad hoc boundary conditions. 
It requires self-consistent calculations of dust dynamics 
combined with radiative transfer calculations and dust 
sublimation. 
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